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Abstract 

~yj-^ ■ A new approach to the analysis of soft and semihard hadron processes is 

Q\ \ suggested. In the frame of the Quark-Gluon String Model the interaction of 

0^ ■ valence quarks and diquarks and sea quarks (antiquarks) of colliding hadrons 

is taken into account. This one is calculated as the exchange by one nonpertu- 
F ^ [ bative gluon, i.e., the cut-off parameter in the gluonic propagator is included. 

|"t ! ' This one allows us to analyze the inclusive hadron spectra in hadron collisions 



at transverse momenta up to 3 — 4 Gev/c. 



1. Introduction 

As it is known the quark-gluon strings model (QGSM) [1-6] based on the 1/N ex- 
^ 1 pansion in QCD [7-9] had a significant success by the description of different hadron 

characteristics. For example, it reproduces a large amount of existing experimental 
! data about hadron production at high energies [1-6]. 

The QGSM has been successfully applied to the production of hadrons containing 
Q_i' light u, d and s quarks. This model reproduces in great detail and in terms of very 

o ■ few parameters a large amount of existing experimental data about hadron produc- 

tion at high energies in hadron-hadron and hadron- nucleus collisions [1-3]. Besides 
QGSM describes adequately the production of charmed mesons at the explored low 
and moderate energies. 

However the characteristics integrated over transverse momentum p t or at av- 
erage transverse momentum < p t > are considered only in the framework of this 
model. So this model is limited usually by the analysis of soft hadron reactions 
or processes at small transverse momenta p t of produced hadrons. There are some 
versions of QGSM taking into account transverse momenta of quarks in the initial 
hadrons |TD|, [TT], |T2| and []l3j] . It allowed to describe inclusive spectra of hadrons 



produced in hN collisions up to pt — 1. — 1.5 Gev/c. But there was a large sensi- 
tivity of these ones to the initial quark distribution over transverse momentum in a 
hadron which has been parametrised usually in an exponential form. On the other 
hand at large p t the ordinary perturbative QCD can be applied to the analysis of 
hard processes. Semihard ones can be explained by QCD models of type as |L4], [L5[] 
etc. However the question arises whether one can apply the QGSM to the analysis 
of semihard processes if to take into account the colour interaction between colliding 
quarks (antiquarks) and diquarks (quarks) before the creation and dacay of qq (or 
q(qq)) string. 
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So the phenomenological Pomeron exchange corresponding to the cylinder graph 
can be understood in the framework of QGSM [1-3] as the exchange of two gluons 
jn| [17]]. However the perturbative calculation of the elastic quark-quark scattering 
amplitude through two-gluon exchange shows a singularity at t = 0. Although 
this singularity can be cancelled when incorporating quarks in the proton wave 
function, this procedure is not able to reproduce the t-dependence of the differential 
cross section observed in experiment Meanwhile some time ago the Pomeron 
was described by the exchange of two nonpertubative gluons JI9L |2(], f?T]], where 
by nonpertubative it means a gluon whose propagator does not show a pole at 
k 2 = 0. In [22j the two-gluon exchange model was applied to the description of 
nucleon-nucleon scattering where a gluon propagator regularized by a dynamically 



generated gluon mass, derived by Cornwall some years ago [23 



The production of hadrons at high energies is described in the QGSM by "cut- 
ting" the forward scattering diagrams of the "cylindrical" type [1-3]. Each cylinder 
corresponds to the exchange of a single Pomeron. Such Pomeron exchange as the one 
of two non-perturbative gluons can be considered for the calculation of cylinder-cut 
graphs in the framework of the QGSM. This paper is dedicated to the exploitation 
of this idea to the analysis of hadron processes in the framework of QGSM. 



2. General formalism 

Consider the inclusive spectrum of hadrons produced in nucleon-nucleon collision 
in the framework of QGSM. As is known the main contribution to such processes 
at large energies gives the graphs of cylinder-cut type in the s channel correspond- 
ing to the Pomeron exchange in t channel [1-6]. The hadron production can be 



considered in the following manner [II]]. Each of two colliding nucleons is divided 
into a quark and a diquark with the opposite transverse momenta. After the colour 
interaction (by means of non-perturbative gluon) between the quark of first proton 
and diquark of another one, and diquark of the first proton and quark of second one 
two quark-gluon strings are created in the chromostatic constant field; then these 
two strings decay into secondary hadrons. This process is repeated n times during 
the production of n Pomeron showers (or 2n qq chains). Such division of the trans- 



verse momentum considered in |Tl[ is analogous to the so called consequent division 
of energy between n Pomeron showers in the multiperipheral model of the hadron 
production 0. 

In principle, there can be so called regular division of energy and transverse 
momentum between 2n qq chains or n Pomeron showers. But at all versions of 
QGSM nobody considered the colour interaction between quark an antiquark or 
diquark and quark respectively. For the soft processes, i.e., at small p t , quark 
distributions over x on the ends of q — q strings are found from Regge asymptotic of 
graphs (cylinder or planar type) at x — > 1 [1-5]. The decay of each this string into 
hadrons is described by fragmentation functions, the form of which is determined by 
the Regge asymptotic of these graphs too [1-5] . The dependence of hadron inclusive 
spectra on p t is the result of the inclusion of internal hadron transverse momenta 



of quarks only [JTU, |TT] and JI3"| . This approach is good to be applied to hadron 
production at p t < 1. - 1.5 GeV/c 0, 0, 0. 

However if we want to analise the hadron production processes at not large pt, 
for example, p t < 3 — 4 Gev/c, the question arises whether one can use the QGSM 
for this aim. In order to do it we must know the distribution of quarks (antiquarks) 
and diquarks over transverse momenta k t on the ends of qq string after some colour 
interaction between constituents of colliding hadrons. For this aim one can use the 
approach suggested in [22 where the Pomeron is described by the exchange of two 
nonperturbative gluons and it means that a gluon propagator has not a pole at 
k 2 



pOfl . For example, in 



[2TI an approximate Shwinger-Dayson [SD] equation 
for the gluon propagator in the axial gauge has been solved and found a solution 
whose infrared behaviour is less singular than a pole at k 2 = 0. This approach has 
given reasonable agreement with experimental data. Recently in [22| it has been 
shown how a gluon propagator regularized by a dynamically generated gluon mass, 
derived in [23 some years ago, successfully describes nucleon-nucleon scattering in 
two-gluon exchange model. The gluon propagator represented as an approximate 
solution of the SD equation and it had not the infrared singularity at k 2 = as a 
result of a some cut-off parameter itlq or so called "effective" gluon mass. In fact 
the model of ref. |22] has a single parameter which was taken to be the ratio of tuq 
and A = A QC d- 

In this paper we take into account the interaction between quark and antiquark 
(or diquark and quark) before the creation and the decay of the corresponding string. 
The main contribution to this one results in the one gluon-exchange between them, 
which can be calculated using the gluon propagator with the cut-off parameter as 
like as in p2| p3| . So this interaction can be considered as the elastic qq or quark- 
diquark (q(qq)) scattering calculated in the one-gluon exchange approximation when 
transverse momenta of quarks or diquarks are changed only. If one represents the 
distribution of quarks or diquarks in a nucleon over x and pt in the factorized form, 
i.e. 

fr(x,h) = fr(x)g?\k t ), (1) 

where r means the flavour of a quark, x and k t are its fraction of the longitudinal mo- 
mentum and the transverse one respectively; f T (x) is the x-distribution and gf\kt) 
is the /^-distribution of this quark (antiquark or diquark) in the initial hadron. Then 
the factorized form of f T (x,k t ) will be true after the one-gluon exchange, but the 
distribution g^\k t ) will be changed only. After such interaction this one can be 
calculated by the following manner: 



gW(k 1>t ) = C X J DVtefSO^&.t - h~ q t )d 2 q t d 2 k t 



(2) 



or 



g?\k u ) = d I D 2 [{k u - k t ) 2 }gW(k t )d 2 k t 



where D(q 2 ) is the regularized gluon propagator and according to [22, it can be 
written in the following form: 



where: 



2/ 2 

m (q 



/n((g 2 + 4m 2 )/A 2 ) 



-12/11 



(4) 



(5) 



ln(4ml/k 2 ) 

Here in (4), (5) b = (33 - 2n / )/487r 2 ; m = 0.37 Gev for A = 0.3 Gev. As in our 
previous paper [II], |T^] we choose the quark ^-distribution in the initial hadron in 
the Gauss form normalized to 1, 



7___/ 



0o(fct) = -exp(->yk 
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(6) 



with 7 to be a free parameter, 7 = 1/ < A; 2 > where < A; 2 > is the average value of 
the squered transverse momentum of quarks in a nucleon. The constant C\ in (2) 
is determined by the normalization condition: 



(7) 



As it is seen from (3) the gluon propagator, D(q 2 ), at very large q 2 returns to the 
usual one asymptotically, i.e. D(q 2 ) oc q~ 2 , because m 2 (q 2 ) — > at q 2 — > 00. This 
method of the calculation of the /^-distribution of a quark (antiquark or a diquark) 
on the ends of qq string allows us to construct inclusive spectra of hadrons produced 
from its decay as the function of x and p t . Note the sensitivity of calculated 
using (2) to the choice of the form is very small if we consider not small k t . It 
is caused by the strong dependence of the gluon propagator (3) on q t . 

Let us return to this construction in the framework of the version of QGSM 
taking into account transverse momenta of quarks [[11], The expression for the 
invariant inclusive hadron spectrum for the reaction pp — > hX corresponding to 
graphs of the cylinder-cut type [1-6] can be written in the following form Jll| : 



p pp ^hx(x,pt) = E 



da a 
d 3 p 



n=l 



(8) 



where a n is the cross section for production of the n pomeron chain (or 2n quark- 
antiquark strings) decaying into hadrons, calculated in the frame of the "eiconal 
approximation" pi[ , (p n (x,p t ) is the x and pt-distributions of hadrons produced in 
the decay of the n pomeron chain. These functions, 4> n (x,p t ), were represented in 



the following form JD]: 



4>n(x,p t ) 



with 



1 ,1 

dxi / dx 2 ^n(x n ,pt;x 1 ,x 2 ); 

x_|_(n) Jx—(n) 



<S> n {x n ,p t] x u x 2 ) = F^\x + (n),p t ;x 1 )F^\x4n),p t ;x 2 )/F qv (n \0,p t ) + 



(9) 



F^\x + (n),p t ;x 1 )F^(x4n),p t ;x 2 )/Fj n \0,p t ) + 



2(n - 1)F£J (x+(n), p t ; Xl )F® (x_ (n),p t ; x 2 



where x± = x±(n) = 0.5(yxf + x^±x), x t = 2y (ml +pf)/s ; is the mass of the 
produced hadron, s is the total energy squared of colliding protons in c.m.s. Note 



that Feynman variable x in n-Pomeron shower depend on n | IT[| : x n = x/(1 — Xq) 
x ~ 0.25. 



n-1 

j 



The functions F^{x±{n) ) p t ) are represented by convolutions 

F^(x ± (n),p t ;x h2 ) = J d 2 k t f^(x lt2 ,k t )G T ^ h (x ± (n)/x h2 ,k t ;p t ), (10) 

F^(0,p t ; Xl , 2 ) = [ d%f^(x l!2 ,k t )G T -> h (0,Pt) = G T ^ h (0,p t ). (11) 



Here r means the flavour of the valence (or sea quark) or diquark, f^(x, kt) is the 
quark distribution function, depending on the longitudinal momentum fraction x 
and the transverse momentum k t in the n-Pomeron chain: 

G T ^ h (z,k t ;p t ) = zD T ^ h (z,k t ;p t ), (12) 

D T ^h{z, k t ;pt) is the fragmentation function of a quark or diquark of flavour r into 
a hadron h. 

The functions f} n \x) and g( n \k t ) can be calculated by using regular |TU 



or con- 



sequential [TT|] division of x and k t between 2n qq chains. Here we use the second one 
|TT|j , so the g^ l \kt) are calculated according to JTTJ] by the following manner. The 



quark functions f T (x, k t ) in initial hadrons are represented in the factorized form (1). 
The distributions functions over transverse momentum p t , g^ikn^)-, are calculated 
analogous to (2). This means that the quark distribution after one nonperturbative 
gluon exchange is given by (2). After exchange by the second gluon the quark dis- 
tribution in corresponding chain will be expressed already via the function g^(k t ), 

#(M = C 2 J g^\k ht )D 2 [{k u - k ht ) 2 ]d 2 k ht . (13) 

Repeating this iteration procedure we obtain the quark distribution function in the n 
chain expressed via the function g^ 1 ' (k n _i >t ) and therefore via the function g^(k t ): 



gi n \k n , t ) = C n J g^\k n ^ t )D 2 [{k n>t - k n ^ t ) 2 }d 2 k n _ ht = (14) 
J d 2 k n _ u D 2 [{k n ^-k n _^ t ) 2 } J d 2 k n _ 2it D 2 [(k n _ ht - k n _ 2it ) 2 ]... 

'd 2 k t D 2 [(k ht -k t ) 2 ]giO)(k t ). 



It is easily to see that in the n-pomeron chain the quark functions will be fac- 
torized too: 

f^(x n ,k n , t ) = fi n \x n )gi n \k n , t ), (15) 

This calculation of ^-distribution of quarks on the ends of 2n-th sring is a principal 
difference from our previous papers [JIT], |T2[. But x-distribution after n-Pomeron 
exchanges f™{x) is calculated according to [PD , i-e. the Feynman variable x in 
n-Pomeron shower depend on n (see above). 



The functions G T ^h{z, kt]Pt) have represented in the following form fTT 

G T -> h (z,k t ;pt) = G T ^ h (z,p t )g T ^ h (k t ), 

where 



1 

9r^h{h) = -exp{-^k\ 

7T 



k t =p t - zk t , 



(16) 

(17) 
(18) 



Substituting now the functions (14)-(17) to (10) and integrate it over d 2 k t . In 
principle, it can be made by analitically if the function g^\k t ) calculated by (10) is 
approximated by a sum of exponantials over k\ (see Appendix). 



3. Results and discussions 

The results of calculation of differential cross sections performed with the help of 
the modified QGSM are presented in Figs. 1-5. Using (8) we can calculate different 
characteristics of hadrons produced in N-N collisions, for example, correlation of the 
average transverse momentum < p t > and hadron multiplicity iV (see for example 
|lT|l). Besides the inclusive spectra of hadrons defined by (8), we can calculate the 
distributions of hadrons over x or pt,i.e.: 

F(x) = J p pp ^ hX (x,pt)d 2 pt (19) 

and 

da ir i- r dx 

^2 = 2 V s J e P P->hx(x,Pt)-j^ (20) 

where E* is the energy of final hadrons in N-N c.m.s. 

The results of the calculations of < pt > (N) for pp collisions are shown in Fig. 1. 
The calculations reproduce the data at small multiplicities N and go below at large 
N. The discrepancy between the results of the calculations and the experimental 
data about < p t > (N) may be connected with the fact that the expressions for the 
cross sections o n (see (8)) were obtained with inclusion of only non-enhanced graphs 
of Reggeon theory [2] . Inclusion of enhanced- type diagrams leads to the appearance 
of 1/x terms in the distributions for sea quarks, whose contribution is especially 
large for x ~ 0. 

The calculation results for inclusive hadron spectra performed using the formulae 
(8), (20) are presented in Figs. 2-5. In Figs. 2,3 we show the inclusive spectra of 
7r + and 7r~ mesons produced in pp collisions at pi a b = 200 GeV/c as a function of 
p t for three values of the Feynman variable x: x = 0, 0.3 and 0.6. We see a good 
description of experimental data with modified QGSM. Corresponding parameter 
values in formulas (6) and (17) are: 7 = 9 (GeV/c)~ 2 what corresponds to < k\ >= 
I/7 ~ 0.1 (GeV/c) 2 , 7 = 7GeV~ 2 . Note that the calculation results depend on the 
value of this parameters very weakly. 



In Fig. 4 we show the invariant cross section E(da/d 3 p versus transverse mo- 
mentum p t for the 7T° for mesons produced in pp collisions at yfs = 52.7 GeV. From 
this figure it can be seen that our version of the QGSM give good description of 
the data [27] up to p t < 3 GeV/c. Some discrepancy with data at larger p t can be 
caused by the following. First of all the p t quark distribution in initial hadrons has 
been taken in the simple Gauss form (6) which is true at small p t . Second reason 
can be conjugated with the Regge trajectories, a(t), which are taken in the QGSM 
at t — 0. We suppose that in the calculations it is needed to take into account the 
dependence a(t) on transfer momentum t. 

In Fig. 5 we present the calculation results for the differential cross sections 
da /dpi of the reactions pp — > D n X, n = + (a), — (b), (c), (<f), at yfs = 27. 4 
GeV. The curves correspond to the case when the Regge trajectory has intercept, 
a>y(0) = 0,, which corresponds to the nonlinear \l/-trajectory. The respective values 
of parameters in fragmentation functions are []12" |: ao — 0.1 10~ 3 , a± = 5. From this 



figure it can be seen that the version of the QGSM under consideration give good 
description of the data [28]. 

The inclusive spectra of all D-mesons have been calculated with the parameter 
value 7, 7 = 2 (Gev/c)~ 2 . The better coincidence of our calculations with the 
experimental data is for D + and D° mesons. 



4. Conclusion 

In this work we have discussed one possible modification of the Quark-Gluon String 
Model to the description of semihard processes. Semihard processes can be explained 
by the QCD models such as [|TJ], [D| etc. The main question considered in this work 
is whether one can apply the QGSM to the analysis of semihard processes if one 
include the colour interaction between colliding quarks (antiquarks) and diquarks 
(quarks) before the creation and dacay of qq (or q(qq)) string. 

Our approach to the analysis of soft and semihard hadron processes in the QGSM 
takes into account the dependence of quark distributions in hadrons and the quark 
fragmentation functions on the transverse momentum p t . In the framework of the 
QGSM the colour interaction of valence quarks and diquarks and sea quarks (an- 
tiquarks) of colliding hadrons has been taken into account. This one has been 
calculated as the exchange by one nonpertubative gluon, i.e., the cut-off parame- 
ter in the gluonic propagator has been included. The method proposed to include 
the p r dependence in the QGSM is close to the successive division of the transverse 
momentum p t and x between n-Pomeron showers or 2n quark-antiquark chains con- 
sidered in our early works [11], [12). This method shows a strong dependence of 



hadron characteristics on n. 

The QGSM modified by such a way has allowed us to analyze the correlation 

< p t > (N) and inclusive hadron spectra produced in hadron collisions at transverse 
momenta up to 3 — 4 Gev/c (Figs. 1-5). The calculations reproduce the data on 

< Pt > (N) at small multiplicities N. To describe the data at large N one need, as 
was mentioned above, to take into account the graphs of enhanced type of Reggeon 



theory [2]. Inclusion of such diagrams result in the appearance of 1/x terms in 
the distributions for sea quarks, which give especially large contribution for x ~ 0. 
Modified QGSM can reproduce a lot of data about invariant cross sections versus 
p t , for example of mesons produced in pp collisions at different values of x. 

In the framework of new version of QGSM, the description of the differential cross 
sections of D + , D~, D° and D° mesons versus transverse momentum p t has been 
performed. The comparison of the calculations with data show good agreement for 
the case when Regge \l/-trajectory has intercept, a^{0) = 0. Such intercept value 
corresponds to the nonlinear \l/ trajectory. Main uncertainties of the calculations 
comes from poor knowledge of intercept of \l/-trajectory. At high energies, the 
predictions for a*(0) differ by a factor close to 3. Unfortunately , the large error in 
charm cross section measurement at = 630 GeV does not allow one to extract a 
useful constant on a^(0). 

The authors would like to thank Yu.A. Budagov for support and constant inter- 
est to this work, A.B.Kaidalov, O.I Piskunova and K.A.Ter-Martirosyan for useful 
discussions. 

5. Appendix 

Show if the function D 2 (q 2 ) can be approximated in form: 

D 2 {q 2 t ) = Y j a t exp{-b t q 2 t ), (21) 
i=i 

where aj,6j-some parameters, then the integration of (10) is performed by analiti- 
cally. Consider firstly the simple case n — 1, then substituting (21) and (6) to (2) 
we have: 

7 N f ^ 
g^\k u ) = - J^ai / exp(-^k lt )exp(-biql)5 {2] \k t - k lt - q t )d 2 k u d 2 q t . (22) 

71 i=i J 

Integrating (22) over d 2 q t and d 2 k u we get: 

N n 

g^\k 2 )=Y.^d l exp{-d l k 2 ) ) (23) 

i=l °i 

where di = Substituting now (1) and (16) to (10), we get the following expres- 
sion for F^(x±, p t ; Xi i2 ) '■ 

F^(x ± ,p t ;x 1:2 ) = fM(x lt2 )G T ^ h (z,p t ) j g^\k t )g T ^ h {p t - zk t )d 2 k t (24) 

Substituting now (23) and (17) to (24) we integrate it over d 2 k t and finally have: 
F^(x ± ,p t ; Xl , 2 ) = f^ix^G^^p^hiz,^), (25) 



where: 

N CLi 

h{z,p t ) = T7zi ex P(-7ziPt), (26) 

here 7^ = d ^J 2 - But for any n (10) is integrated over d 2 k t anologous to this simple 
case. 

The question arises whether D 2 (q 2 ) (see (3)) can be represented in the form (18). 
This approximation can be true at the following parameters: N = 5, a\ = 0.31, a 2 = 
0.333, a 3 = 0.172, a 4 = 0.155, a 5 = 0.026; and b x = 9.294, b 2 = 2.78, b 3 = 3.619, b 4 = 
0.882, 6 5 = 184. 
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FIGURE CAPTIONS 
Fig. 1. Dependence of < pt > on the number of charged particles, N c h- The 
curve shows our calculation for yfs = 63 GeV. Symbols denote the experimental 
data for y/s = 63 GeV, rapidity interval, \y\ — 2 



Fig. 2. Invariant cross section of ir + mesons produced in reaction pp — > n + X at 
Pm — 200 GeV/c as a function of p t for the values x = 0, 0.3 and 0.6. Experimental 
data are from Refs. [^6 . 



Fig. 3. Invariant cross section of tt~ mesons produced in reaction pp — » tt^X at 
Piab = 200 GeV/c as a function of p t for the values x = 0, 0.3 and 0.6. Experimental 
data are from Refs. [26|. 

Fig. 4. The inclusive cross section of the reaction pp — > tc°X versus transverse 
momentum at ^/i = 52.7 GeV. The data are from Ref.[|27|. 

Fig. 5. Differential cross sections of the reactions pp — > D + X (a), pp — > D~X 
(b), pp — > D°X (c), pp — > D°X (d) at the energy -^i = 27.4 GeV as a function of 
the square p\ of the transverse momentum. Solid curves correspond to calculations 
with the intercept a^,(0) = 0. The data are from Ref. 



